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TECHNICAL NOTE 3832 


ON POSSIBLE SIMILARITY SOLUTIONS FOR THREE-DIMENSIONAL 
INCOMPRESSIBLE LAMINAR BOUNDARY LAYERS 
II - SIMILARITY WITH RESPECT TO STATIONARY POLAR COORDINATES 
By Howard Z. Herzig and Arthur G. Hansen 


SUMMARY 

Solutions of mainstream flow patterns for three-dimensional,, laminar, 
incompressible thin-boundary- layer flows (over flat or slightly curved 
surfaces) having similarity with respect to stationary polar coordinates 
in the plane of the surface are derived. The solutions are summarized in 
a table. 


INTRODUCTION 

Considerable attention has been devoted in laminar-boundary- layer 
research to theoretical solutions of the two- and three-dimensional 
incompressible-boundary- layer equations lasing the "similarity" approach. 

In this method, the partial differential boundary- layer equations are 
transformed by means of a similarity parameter Tj and rewritten In terms 
of functions of T), their derivatives, the mainstream velocity components, 
and their derivatives. Solutions are then sought for the mainstream flow 
conditions for which the transformed equations reduce to ordinary differ- 
ential equations for the functions of T} (refs. 1 to 10 ). Some experi- 
mental evidence is presented In reference 10 in support of this kind of 
theoretical development for laminar flows. Reference 11 presents a sys- 
tematic approach to similarity- type solutions using a generalized simi- 
larity parameter. As a result, reference 11 has obtained solutions for 
the permissible mainstream flows for all the boundary- layer flows having 
classical similarity with respect to stationary rectangular coordinates. 
The present report is an extension of the work of reference 11. Solutions 
are sought for the mainstream flows in stationary cylindrical coordinates 
for all the boundary- layer flows having classical similarity with respect 
to the polar coordinates in the plane of the surface. 
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SYMBOLS 


a,b, c,C 

F, F(n) 

f 

G, G(ti) 
g>g(r,e) 
k,m,n 

r,Q,J 

U,W 


u.,v,w 


W 

n 

v 

Subscripts: 


constants ~ 

function of similarity parameter, u = UP" (t)) 
arbitrary function 

function of similarity parameter, w = W3*(n) for W / 0 

function of coordinates r and 0 

constants 

polar coordinates 

mainstream velocity components in 0, r directions, 
respectively 

boundary- layer velocity components in 0,y,r directions, 
respectively 

function of coordinates r and 0, w s WG’(ti) for W = 0 
similarity variable, rj = yg(r, 0)//\/v 
coefficient of kinematic viscosity 


1,2,3, . . . index numbers 

Primes denote differentiation. 


ANALYSIS 

Boundary- Layer Equations in Stationary Polar Coordinates 

The three-dimensional laminar, incompressible, thin-boundary- layer 
equations in stationary cylindrical coordinate form for flows over flat 
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(or slightly curved) surfaces with coordinate axes as shown in this sketch 
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are given by 


uw 

u 3u 

r 

+ 7 Se 

u 2 

u dw 

r 

+ x'Se 


c)u . 

<3u 

S 2 U 

W ^ — + v ^ — — 

or oy 

V Sy2 

dw 

3w 

S 2 w 

w ^ — + V ^ — — 

or oy 

V By 2 


UW U dU __ du 

~ + r ^0 W 3r 

U 2 U Bw TT dw 


(la) 


(lb) 


where u, w, and v are the boundary- layer velocity components in the 
0, r , and y directions, respectively. 


Equations (la) and (lb) are the boundary- layer flow equations in the 
tangential and radial directions, respectively. Consistent with the 
restriction to thin-boundary- layer flows for flat (or slightly curved) 
surfaces, the mainstream velocity components are 


U = U(r,0) 
W = W(r,0) 
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The continuity equation for the boundary- layer flow is given by 

1 du bv w , c*v . _ \ 

?Se + 5? + ? + 5F- 0 < lc) 

The boundary conditions are ~r.. 

u = w = v = 0 for y = 0 

u -> U^i \ (Id) 


as y 


w -+■ ¥ 


Transformation of Equations Using Generalized r| 

The method of search for symmetric solutions described in reference 
11 suggests exact solutions (i.e., solutions for flows, having similarity 
with respect to the polar coordinates r , 0 in the plane of the surface) 
of equations (l) can be obtained as a result - of transformations to new 
coordinates r, 0, and t] where the space variable r[ is of the form 

n 2 g(r,e) (2) 

/ yv 

In the rectangular coordinate systems, when either component of 
mainstream flow equals zero, straight mainstream flows result, and the 
boundary- layer flow is two-dimensional with no secondary- flow overturning. 
In polar coordinates when U = 0, the mainstream flows are straight and 
u = 0 (ref. 12). However, when U ^ 0, there is curvature of the flow 
streamlines even though ¥ = 0, and three-dimensional boundary- layer 
overturning does result. In polar coordinates, as will be seen subse- 
quently, it will be necessary to treat separately the cases of ¥ = 0 
and ¥ / 0. 

¥ / 0. - Following the rectangular- coordinate- system analyses 
(ref. 11 ), 


u = UF’ (tj), U/O 
w = WG'(n)> ¥ / 0 


(3a) 

(3b) 
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The conditions on F* and G' required to satisfy boundary conditions 
on u and w are 


F*(0) = G'(0) = 0 ^ 
lim F'Cq) = 1 


> 


lim G’Cn) = 1 

f] ->• °o J 


(4a) 


Now v may be determined by integration of the continuity equation 
using equations (3): 


v = 


_ - 'V/v 


g 


(l Su U S In g \ fbv TT a In g . w\ „ 

\r - r / F + - W §r + r j G 




( 5 ) 


In order that v = 0 for y = 0 as required, it is possible with- 
out loss of generality to set the boundary conditions 


and 


F(0) = G(0) = 0 


f(r,e) = 0 


(4=b) 


(See appendix C of ref. 11 for a discussion of the necessary and suffi- 
cient boundary conditions . ) 

Upon substitution of equations (3) and (5), equation (la) becomes 

(W f 0) 


W 


i Su 


S in U 


^ (F‘G* - GF" - 1) + ^ ^ (F' 2 - FF" - 1) + W u - u (F‘G r - l) + 

■gj " - gt — FF " + I ^ ^ GF" - GF" - g 2 F"’ = 0 (6) 


and equation (lb) becomes 


t (1 - *' 2 ) + ? (r .o. - 1) + I? (0-2 - GG" - 1) - i fg FG" 


Wr 


3? 




U 5 lng 2 JQI, + W 5 ln _ g 2 Q(j" - — GG" - g 2 G’" = 0 


lr 


" 5 ? 


( 7 ) 
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The purpose of this investigation is to determine mainstream flow 
solutions for which the transformed equations (6) and (7) reduce to 
ordinary differential equations. As in reference 11, the mainstream flow 
conditions that make the coefficient of the functions of t} proportional 
are sought. Under these ordinary-differential- equation conditions 
(abbreviated to o.d.e. conditions), the common variable terms in the 
equations may be divided out, leaving ordinary differential equations for 
F and G. (The actual numerical solutions of the ordinary differential 
equations are not attempted herein.) Although coefficients of similar 
functions may be grouped and made proportional, the two techniques can be 
shown to be equivalent. 

For convenience, the coefficients for the functions of T) in equa- 
tions (6) and (7) are presented here in the order of their appearance. 
With W f 0, they are: 

© 

© 

© 

<D 

© 


w 

r 

© 

aw 

3F 


1 dU 
r 3© 

© 

g 2 


a In U 
W dr 

© 

u 2 

Wr 


U a In g 2 

© 

u a_ 

In W 

2r 

r 

60 

W a In g 2 





2 ~5F 


The o.d.e. conditions require these nine coefficients to be pro- 
portional to each other. The most general case, W = W(r,0), can be 
solved readily. From o.d.e. conditions on (I), (j) , and (§) , 

W = c-jU (8) 


g 2 = 


c 2 U 


Then from o.d.e. conditions on d> and ® , using equation (8), 


(9) 


dU = f_3 
3© c^_ 


U 


( 10 ) 


which may be integrated to give 


U = 


f5e 

e 1 ^i(r) 


( 11 ) 


4177 



4177 


NA.CA TN 3832 


7 


From o.d.e. conditions on © and (§) , using equations (8) and (ll), 

q(r) 


1 

= c„ — 


fi (r) “ u 4 r 


( 12 ) 


Therefore, 


fl(r) = c 5 r 


(13) 


Using equation (9) 
permits W, U, and 


and redefining the constants involved for convenience 
g 2 • to be written 


W(r, 0 ) = br n e mS 
U(r,9) = ar n e me 
g 2 (r,S) = cr n-1 e m ® 


(14) 


O 

No further restrictions on the forms of W, U, and g^ are required 
to satisfy the o.d.e. conditions for the remaining coefficients 
© , and Substitution of equations (8) and (9) indicates that © 

and © are already proportional to © , while © and © are proportional 
to © . 


In addition, further analysis shows that the remaining possible main 
flows W = W(r), W = W(0), or W = constant / 0 can be obtained from 
equations (14) by suitable choices of n, m, a, and b. For example, 
when m = 0, 



(15) 


The ordinary differential equations for W(r,0) ^ 0 now can be 
obtained by substitution of equations (14) Into equations (6) and (7)j 


b(n + l)(F‘G' - l) + am ( f ' 2 - - 1^ - 


GF" - 


GF - cF m = 0 


(16) 


~ (1 - F' 2 ) + am 


^F‘G‘ - - 1^ + bn (G* 2 - 1) - GG" - cG m = 


( 17 ) 
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W = 0 . - When W = 0, the corresponding boundary^- layer equations 
(l) become 


uw u du 

du 5u 

: d 2 u 

U du 

(18a) 


V dy 2 = r-35 

o 

u ^ u bw 

bw bw 

' d 2 w 

u 2 

(18b) 

■~ + r^e +v ^ + v ^ 

dy 2 


The equation of continuity for the boundary- layer flow remains unchanged; 

(lc) 


1 5u , Sw , w , Sv - 
r3e + 3? + r + ^ = 0 


The boundary conditions now are - — • 

u = w = v = 0 for y = 0 | 

u -> > (18c) 

V as y oo | 

w -*■ 0 J J 

For main flows such that W = 0, u and w are defined as follows: 

u=UF’(ti) (19a) 

w=WG'(ti) (19b) 

where W = W(r } 8) ^ 0. The boundary conditions on FV_ and G f required 
to satisfy boundary conditions on u and w in equations (18) 

(W = 0, W(r,0) / 0) are 


F'(0) = G'(0) = 0 
lim F* (rj) = 1 

Tj -*• * 




> 


lim G* (tj) = 0 

TJ -*• oo J 


(20a) 


The expression for v obtained by integration of the continuity 
equation (lc) is the same as equation (5) with W being replaced by Tf; 


- - 1 *+(§ - - 
yp- - yG' + f(r,e) 


( 21 ) 
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As "before, the "boundary conditions chosen as sufficient to provide 
that v = 0 for y = 0 are 


F(0) = G(0) = 0 

f(r,e) = 0 


(20b) 


Substitution of equations (19) and (21) into equations (18a) and (18b) 
produces 


? (F*G* - OF" ) + i ^ (F* 2 - FF" - 1) + W 5 - ^ - U F‘G« + FF" + 

2 — GF" - GF" - g 2 F m = 0 


r 3e 

W 5 In g 2 


5w 


( 22 ) 


t- 

Wr 


(1 - F l2 ) + 2 a f*G' 


Sw 

37 


(G 


‘ 2 " GG ”) “ 7 + 


■gj 5 FG" + | S g - GG" - | GG" - g 2 G'" = 0 (23) 


The argument concerning determination of the o.d.e. conditions by 
means of relations between the coefficients of the functions of rj in 
equations (22) and {23) remains unchanged. These coefficients, it may be 
noted, are the same as the coefficients for equations (6) and (7) with W 

replaced by W- Thus^ the expressions for W, U, and g can be -written 

W(r,e) = br n e m8 for b f 0 ^ 

W = 0 

U(r, 8) = ar n e m0 
g 2 (r,e) = cr n-1 e m0 J 

The resulting ordinary differential equations are (W = 0) 
b(n + 1) F’G* + am ^F t2 - ^ - 1^ - * - 5 ^ GF" - cF'" = 0 (25) 




(24) 


(l - F' 2 ) + am 


^F'G* - + bnG' 2 - - 3 ) . GG" - cG m = 0 


( 26 ) 
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RESULTS AND DISCUSSION 

The analysis of three-dimensional, laminar, incompressible, thin- 
boundary-layer flows having similarity with respect to polar coordinates 
has led to solutions for mainstream flows described by equation (14) or 
(24). As a result of this analysis, table -I has been - prepared, which 
summarizes the two cases of mainstream flows over a flat or slightly 
curved surface for which the boundary-layer flows have the required 
similarity. 

As described earlier, secondary flows exist even when the radial 
.component of mainstream flow vanishes; W = 0. For such cases a function 
W = W(r ,0) ^ 0 is defined and the boundary- layer radial component of 
flow is expressed as 


w = WG»(t]) (19b) 

The table presents these two cases W / 0 and W = 0, W f 0, for the 
corresponding forms of the tangential component U = U(r,6), which permits 
a solution by reduction of the boundary- layer equations to ordinary dif- 
ferential equations. 


Mainstream 

When W f 0, the mainstreams are spiral flowB. For W = 0, circular 
mainstream flows are obtained. 

U,W, W . - In regions where the thin-boundary- layer theory is applica- 
ble, the mainstream is very nearly parallel, to the surface; U and W 
are functions of r and 0 only. 

The analysis shows that for_the similarity solutions considered here 
only one form of U and W (or W) ie possible; that Is, 

U = ar n e m0 

W (or. W) = br n e m ® 

When W = 0, by choosing m = 0 (case I), many of the most frequently 
encountered flow distributions can be obtained by suitable choices of n. 
Corresponding to n = -1 is free-vortex flow; for n == 1, wheel-type 
flow. In all cases, W/U or w/u are constants. 
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Projection of mainstream on surface . - The equation for the projec- 
tion of the mainstream on the surface may he obtained by integrating 


W _ dr 
U “ r d0 


(27a) 


for each case yielding 


r 


£ e 


c c e 


(27b) 


The slope of the projected streamline with respect to 0=0 obtained from 

r + 4s tan S — tan 0 + 1 

slope = , — = (28) 

-=-q - r tan 0 — - tan 0 __ _ 

d0 a 

is found to be independent of radial position r. 

Irrotationality. - For the mainstream flows considered in this in- 
vestigation and in regions of thin boundary layers (as assumed for the 
analysis), only the component of vorticity normal to the surface 

1 Sw SU U 

r 30 ' SF ‘ r 

can be much different from zero (ref. 11 ). The values of the constants 
specified under the chart listing "Irrotationality" were obtained i n each 
case from 

i - |2 - f - (mb - a n - a) = 0 (29) 

These values serve to set the conditions for nearly irrotational main- 
stream flows. 


Boundary Layer 

As discussed in reference 11, the physical interpretation of the 
boundary- layer behavior that the mathematical representations purport to 
describe is best found by examining the behavior of T] and, in particu- 
lar, g(r,0). The boundary- layer thickness was shown (ref. ll) to be 
inversely proportional to g. In order for the theoretical boundary layer 
to have a beginning at a leading edge, as in a real fluid, there should, 
therefore, be a line along the surface for which g{r,0) is infinite. 
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In the solutions presented here, this occurs in the finite part of the 
plane only at the point r = 0 for values of n < 1. For n > 1, the 
boundary layer may be considered to have a beginning only at r = 
there the mainstream velocities take on "infinite” values. 


Ordinary Differential Equations - 

The actual numerical solutions of the ordinary differential equa- 
tions are beyond the scope of the present investigation. The literature 
contains examples of numerical solutions that have been calculated for 
particular values -of the constants. Some of these examples are noted in 
the listing "Comments and references" associated with each case in the 
chart . 

The present analysis then serves only to display the ordinary dif- 
ferential equations that can be obtained with the underlying assumptions. 
In any particular case of interest for which the equations are appropri- 
ate, the existence of the numerical solution and its computation must be 
obtained individually . Nevertheless, some general remarks can be made 
here (as in ref. ll) concerning the numerical solutions. 

Separation of F and G . - Under certain choices_of the free con- 
stants involved, the functions F and 0 are separable; that is, one 
equation of the pair of ordinary differential equations will contain 
terms in only one of these functions and its derivatives. Numerical 
solutions are much more readily obtained in such cases than when the 
functions are not separated. It can be noted from the table that by 
choosing a = 0 in case I, equation (17) contains only terms in G and 
its derivatives. This corresponds to mainstreams having no tangential 
components of flow. Equation (17) becomes a Falkner-Skan type equation 
with known solutions (refs. 1 and 2). Although it is not apparent from 
the equations in the table alone, when a = 0, in case I, then u = 0 
(ref. 12), equation (la) disappears, and bo ioes equation ( 16 ) . Such 
flows are really two-dimensional flows originating from a stagnation 
point and flowing out along straight radial lines. 

Linearity in u or w . - As discussed in reference 11 and applied 
in reference 10, an extension of the solutions beyond strict similarity 
of the velocity component can sometimes be made by addition of solutions 
where the boundary- layer equations are linear in u or w. Such exten- 
sions are not possible for the boundary- layer flows investigated here, 
because equation (lb) is always nonlinear in u and in w except for 
the trivial case of no mainstream flow. — 
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CONCLUDING REMARKS 

Solutions are obtained for the mainstream flow patterns for boundary- 
layer flows having classical similarity with respect to stationary polar 
coordinates. The results are summarized in the table. The exact solu- 
tions obtained are beset with the difficulty that their boundary layers 
have no proper leading edge in the finite part of the plane, whereas in 
turbomachines a definite leading edge is generally required. Neverthe- 
less, the analysis enables a study of the properties of the boundary- 
layer flows and may have direct applicability when attention is confined 
to appropriate regions of the flow. 

The solutions are considered completed for the sake of this investi- 
gation when the boundary- layer equations have been transformed into ordi- 
nary differential equations. The actual numerical solutions for the 
ordinary differential equations so derived are beyond the scope of this 
work . 


There are three configurations of main- flow streamlines for which 
the similarity solutions here could be obtained. The mainstream may be 
(l) a stagnation- type flow out along radial lines from a stagnation 
point, (2) spiral flow out from (or in toward) a central point, or (3) 
circular flow. By suitable choice of the free constants involved, the 
flows may include at will cases of acceleration or deceleration in the 
radial or tangential directions. The solutions obtained here thus con- 
stitute an extension of the similarity solutions obtained in reference 11 
for rectangular coordinate systems. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 24, 1956 
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TABLE I. - SIMILARITY SOLUTIONS IN STATIONARY POLAR COORDINATES 

CASE I 


V / 0 


u 

ar n e m0 

w 

br n e™*, b / 0 

Tt 

/ or n-l e »e\ l/2 / c0 \I/2 / cvV/2 , 0 

y\ y ) ylyva r) - y^J * c r 0 

Ordinary 

differential 

equations 

(16) b(n + 1)(F'G' - 1) + am(p' 2 - ^ - lj - OF" - oF"' - 0 

(17) -^(1 - F' 2 ) + am^F'Q 1 - %jr - l) + bn(0' 2 - l) - 30" - oG m - 0 

Boundary 

conditions 

P'(O) - Q ' (0) - P(0) - 0(0) - 0; lim P'(D) ■> llm 0'(T?) = 1 

7?-t« 7J-t« 

Projection of 
mainstream 
on surface 

| 0 

r » Ce , spiral flow streamlines, (a / 0) 

Irrotation- 

ality 

bin - an - a « 0 

Linearity 
in u 
Linearity 
In w 

(la) ra = 0 

(la) Linear 

( lb ) Always nonlinear 

Separation 
of F and 
G 

(16) Not possible for this case, (b / 0) 

(17) a - 0 

Comments and 
references 

a « 0, stagnation flow. Eq. (16) vanishes. Eq. (17) becomes a Falkner- 
Skan equation, which la completely solved In refs. 1 and 2. 

Ref. 3: a » m - 0, b - o/2, n *■ 1. 

Ref. 13, p. 71: plane stagnation flow: a « m » 0, n •« 3, b « 1/3, 

c *« l, 0 » $ (ref. 13) . 

Ref. 13, p. 74: Three-dimensional stagnation point flow, axisymmetrlcal 

case: a - a ■* 0, n-1, b-c, 0 ■ f (ref. 13). 
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TABLE I. - Concluded. SIMILARITY SOLUTIONS IN STATIONARY POLAR COORDINATES 

CASE II 


W = 0, CT / 0 


u 

ar n e m ® 

w 

0 

T) 

/• cr n-l eraS \ l/2 ( oU ^l/2 

yl v ) y^var/ » c r o 

Ordinary- 

differential 

equations 

(25) b(n + 1 ) (F'G 1 ) + am^F' 2 - - 1^ - b ^ g + GF" - cF m - 0 

(26) 2§-(l - F' 2 ) + am^F'G' - + bn(o' 2 ) - b (?g + .$1 GO 11 - cG m - 0 

Boundary 

conditions 

F ' (0) - G'(0) - F(0) - G(0) - 0; lira P' (J?) - 1; lira G'(TJ) - 0 

T) -*■<*> 

Projection of 
mainstream 
on surface 

v <* C, circular flow streamlines 

Irrotati on- 
allty 

a(n + l) =■ 0 

Linearity 
in u 
Lineari ty 
in w 

(la) m = 0 

(la) Linear 

(lb) Always nonlinear 

Separation 
of F and 
G 

(25) Not possible for this case ~ 

(26) a » 0 

Comments and 
references 

a =■ 0, no flow 

Ref. 13, p. 157, eq. (10.9), (discussing work of ref. 14): 

a ■= b =* e, m =» 0, n-1. F' here = G (ref. 13). O' here - F 
(ref. 13) 


NACA - Langley Field, Va. 
















